
import sys
sys.path.append( '../lib/' )

import numpy as np

import pylab

from gtd import *


pi = np.pi
hw = 7.0
h0 = 4.0
R0 = 13.0

puntos = []

Rr = np.linspace( 1, 45.0, 1000 )

for R in Rr:
	
	phi  = 2*pi - np.arctan( R0/(hw-h0) )
	
	phip = np.arctan( (R-R0)/hw )
	
	s    = np.sqrt( R0**2 + (hw-h0)**2 )
	
	x    = np.sqrt( (R-R0)**2 + hw**2 )
	
	k    = 2*pi/( 0.002 )
	
	puntos.append( np.abs( EdiffPlane( 1.0, s, k, phi, phip ) ) )
	
pylab.plot( Rr, puntos )
pylab.show()
